# -*- coding: utf-8 -*-

R = 8.314

T0=298
T=293

def cp1(a,b,c,d):
    cp = a + b/2*(T0+T) + c/3*(T0*T0+T0*T+T*T) + d/4*(T0+T)*(T0*T0+T*T)
    return cp

def cp2(a,b,c):
    cp = a + b/2*(T0+T) + c/(T0*T)
    return cp

# для теплоемкостей
# этиловый спирт
a=19.07;	b=105.2/1000;	c=-31.04/1000000;	d=0;
print cp1(a,b,c,d), cp1(a,b,c,d)-R

# N2
a=27.87;     b=4.27/1000;	c=0;
print cp2(a,b,c), cp2(a,b,c)-R

# O2
a=31.46;     b=3.39/1000;	c=-3.77*100000;
print cp2(a,b,c), cp2(a,b,c)-R

print '-----'
T=2300

# CO2
a=44.14;     b=9.04/1000;	c=-8.53*100000;
print cp2(a,b,c), cp2(a,b,c)-R

# вода
a=30.0;     b=10.71/1000;	c=0.33*100000;
print cp2(a,b,c), cp2(a,b,c)-R

# N2
a=27.87;     b=4.27/1000;	c=0;
print cp2(a,b,c), cp2(a,b,c)-R